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We present the prospects of searches for neutral, long-lived particles which decay to photons 
using their time of arrival measured with a newly installed timing system on the electromagnetic 
calorimeter (EMTiming) of the Collider Detector at Fermilab (CDF). A Monte Carlo simulation 
shows that EMTiming can provide separation between decay photons from particles with both a 
long lifetime and a low boost, and prompt photons from Standard Model backgrounds. Using a 
gauge mediated supersymmetry breaking (GMSB) Xi ~ * lG model we estimate a quasi-model- 
independent sensitivity using only direct neutralino pair production, and also estimate the expected 
95% confidence level exclusion regions for all superpartner production as a function of the neutralino 
mass and lifetime. We find that a combination of single photon and diphoton analyses should allow 
the Tevatron in run II to easily extend the exclusion regions from LEP II at high neutralino masses 
and lifetimes, and cover much, if not all, of the theoretically favored < 1 keV/c 2 parameter 
space for neutralino masses less than 150 GeV/c 2 . 



I. INTRODUCTION 

The electromagnetic (EM) calorimeter at the Col- 
lider Detector at Fermilab (CDF) Q has recently been 
equipped with a new nanosecond-resolution timing sys- 
tem, EMTiming Q, to measure the arrival time of energy 
deposited (e.g. from photons). While it was initially 
designed to reject cosmics and accelerator backgrounds 
0, we investigate the possibility of using it to search 
for neutral particles 0] with a lifetime of the order of a 
nanosecond which decay in flight to photons. An exam- 
ple of a theory which would produce these particles is the 
gauge mediated supersymmetry breaking (GMSB) model 
Q with a neutralino, Xit as the next-to- lightest super- 
symmetric particle (NLSP) and a light gravitino, G, as 
the LSP. In this scenario the neutralino decays prefer- 
ably (~100%) as x? - * lG with a macroscopic lifetime 
for much of the GMSB parameter space. 

We begin with a study of the properties of events 
where timing can be used to separate between decay 
photons from long-lived particles, and photons produced 
promptly at the collision. A suitable variable to describe 
this distinction is the measured difference between the 
time after which the photon arrives at the face of the 
detector, and the time a prompt photon would virtually 
need to reach the same final position. This time differ- 
ence for a prompt photon, from Standard Model (SM) 
sources, is exactly but is always greater than for pho- 
tons from delayed decays, as in GMSB/Supersymmetry, 
if we neglect the measurement resolution. We call this 
difference As: 
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where tf — U is the time between the collision and the 
arrival time of the photon at the face of the detector, 
and \xf — Xi\ is the distance between the final position 
of the photon and the collision point. The situation is 
visualized in Fig. All four variables can be measured 
by the CDF detector Q and give a system resolution of 

^EMTiming ~ 1.0 IIS @]. 

An important note in our analysis is that photons from 
long-lived particles will usually not arrive at the face of 
the EM calorimeter at the usual 90 degree incident angle. 
This could have serious implications for photon identifi- 
cation. For the purposes of this study we assume that 
this issue can be addressed without significant changes to 
the identification efficiency. We further assume that the 
additional handles such as EMTiming and timing in the 
hadronic calorimeters provide the necessary robustness 
needed to convince ourselves that photons which might 
not pass ordinary selection requirements are indeed from 
our signal source as opposed to sources which could pro- 
duce fake photons and missing transverse energy, like 
cosmics. 

We estimate our sensitivity to two different types of 
new particle production using GMSB models. As a quasi- 
model-independent sensitivity estimate to generic long- 
lived particles we simulate direct neutralino pair produc- 
tion and decay, and examine the dependency as a func- 
tion of both neutralino mass, m^, and lifetime, 75. For 
a "full" GMSB model sensitivity, which means including 
all relevant GMSB subprocesses such that the neutrali- 
nos are part of cascades from gauginos and squarks, we 
allow all SUSY particle production and decay, and again 
vary the mass and lifetime variables. To choose analy- 
sis final states for both we consider three issues: 1) with 
neutralino lifetimes longer than a nanosecond it is possi- 
ble that one or both of the neutralinos leave the detector 
before they decay, 2) with gravitinos or the neutralino 
leaving the detector 1£t should also help separate signal 
from SM backgrounds and 3) to ensure that our predic- 
tions are as reliable as possible we want to use the data 
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Figure 1: A schematic diagram of a long-lived neutralino de- 
caying to a photon and a gravitino in the CDF detector. The 
neutralino emanates from the collision at (xi,ti) and after a 
time t it decays. While the gravitino leaves the detector the 
photon travels to the detector wall and deposits energy in the 
EM calorimeter where its final location x / and arrival time t / 
can be measured. A prompt photon would travel directly from 
Xi to x f which can both be measured. The difference between 
the actual time the neutralino/photon needs, At = tf—U, and 
the time a prompt photon would need, ifi—ii!^ ; s defined as 
As. The SM typically produces prompt photons which have 
As = ns, whereas photons from delayed decays from SUSY 
have As > ns, assuming a perfect measurement. 



selection requirements and background predictions from 
previously published papers by CDF Q and D0 Q. In 
the 1992-1995 collider run (run I) of the Tevatron three 
types of analyses match these criteria: CDF and D0 re- 
sults in 77 + $t E3| exclusive 7 + $t (7 + ffir+ jets) 
from CDF and 7 + $ T + jets from D0 Since 
there are no jets at the parton level in direct neutralino 
pair production, in this case we consider analyses with 
final states 77 + $t and 7 + ftr + jets. To estimate 
the sensitivity for full GMSB neutralino production we 
consider both 77 + Ifir and 7 + I$t + jets analyses. 



For both direct neutralino pair production and full 
GMSB production we quantify the sensitivity for 2 fb^ 1 
in run II using the expected 95% confidence level (C.L.) 
cross section upper limits. Results for both with and 
without the EMTiming system, using kinematics cuts 
only, illustrate the contribution to the final sensitiv- 
ity from kinematic and timing information considera- 
tions ^3] • Finally, we compare the final mass and lifetime 
exclusion regions for a GMSB scenario to direct and indi- 
rect searches from the ALEPH experiment ^3] at LEP II 
and the theoretically favored parameter space from cos- 
mological model restrictions of < 1 keV/c 2 (T3 |. 



II. KINEMATIC PROPERTIES OF EVENTS 
WITH LONG-LIVED PARTICLES WHICH 
DECAY TO PHOTONS 

While the final sensitivity studies use both a full 
physics generation and a detector simulation of the ge- 
ometry and timing resolution, we begin with a study of 
the kinematic properties of events which yield large As 
measurements using a "toy Monte Carlo." For now the 
CDF detector is assumed to be a cylinder, with length 
3.5 m and radius 1.7 m, instrumented with time and po- 
sition detectors of perfect resolution for both the collision 
point and where the photon hits the face of the detector. 
Neutral particles, which we will refer to as neutralinos, 
are simulated as emanating isotropically from the center 
of the detector and emit a photon isotropically after a 
lifetime in their rest frame. For a promptly decaying 
photon the minimum time corresponds to the nearest dis- 
tance to the detector face: ^77^ = 5.6 ns; the maximum 

to the largest distance ^(Um) 2 + (^|a)2/ c = 8.1 ns . 

For pedagogic reasons, neutralinos are simulated with a 
flat momentum and lifetime distribution, i.e. indepen- 
dently any lifetime and momentum have equal probabil- 
ity. We note that to be conservative here, as later in the 
paper, only those neutralinos that decay before the face 
of the detector are considered to have produced a photon. 

Figure shows the measured As versus the event life- 
time of the Xi m the ia b frame, r ev t,L- For As > 10 ns 
there is a roughly linear relation between As and r ev t,L- 
For a fixed T GVt ,L the maximum As (upper bound) occurs 
when the neutralino travels to the farthest corner of the 
detector and then emits a photon backward to the op- 
posite corner. Analogously we get a minimum As (lower 
bound) if the neutralino travels with high momentum 
to the nearest part of the detector and emits a photon 
forward. The latter would look like a usual prompt pho- 
ton event except for the difference in velocity between the 
neutralino and the photon. If the event lifetime is greater 
than the maximum time a prompt photon would need to 
travel to the detector then As is restricted from below 
and As > ns (given that the neutralino decays inside 
the detector). Thus, the spread mainly comes from de- 
tector geometry but with the neutralino momentum also 
contributing to the width. 

Figure shows As versus the neutralino boost for the 
lifetime slice 8.5 ns < r cvti L < 9.0 ns. A low boost (be- 
tween 1.0 and 1.5) allows large As since neutralinos can 
have a larger lifetime without leaving the detector. Neu- 
tralinos with high boost are more likely to leave the detec- 
tor, and even if they do not and their photon is detected, 
it has low As (0 ns < As < 2 ns). Thus, events with 
large As are produced by neutralinos with long lifetimes 
and low boost. 

Next we consider the efficiency for neutralinos to re- 
main in the detector and/or produce a photon with large 
As. Figure0]shows the efficiency, the fraction of all gen- 
erated events that produce photons which pass a given 
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Figure 2: The As distribution as a function of the event life- 
time in the neutralino lab frame for a toy Monte Carlo sim- 
ulation. In general, As is proportional to r cv t,L. At large 
r ov t,L most of the neutralinos leave the detector and are not 
shown here. The spread perpendicular to As ~ r cvt ,L origi- 
nates in variations of the neutralino momentum as well as in 
variations in the travel time of the photon due to detector ge- 
ometry. Essentially, events with large As require a neutralino 
with a long lifetime. 



Figure 3: The As distribution as a function of the boost of 
the neutralino for a lifetime "slice" of 8.5 ns < r ev t,L < 9.0 ns. 
In the region 1.0 < boost < 1.5 neutralinos remain in the 
detector and can produce a large As. Neutralinos with high 
boost, that is high pr, are more likely to leave the detector or, 
if they don't, produce low As. Thus, events with the largest 
As are produced by neutralinos with large lifetimes and low 
boosts. 



III. NEUTRALINO PAIR PRODUCTION AS A 
MEASURE OF QUASI-MODEL-INDEPENDENT 
SENSITIVITY 

A. Analysis Methods and their Efficiency as a 
Function of Neutralino Mass and Lifetime 



As restriction, as a function of the event lifetime, r ev t, 
for As > (neutralino stays in the detector), 3 ns and 
5 ns for the same production distribution. While these 
results change for a more realistic pr spectrum, the qual- 
itative features are instructive. In the limit of T ovt = ns 
and As > the efficiency is 100% and the efficiency de- 
creases with higher event lifetime, since the neutralinos 
are more likely to leave the detector. When one applies a 
As cut however, there is no efficiency for events that con- 
tain neutralinos with a low event lifetime (r ev t < 2 ns). 
For any As > requirement the efficiency goes to 0% 
at Tovt = ns, since all photons would have As = 0. 
A higher As cut gradually suppresses events with a neu- 
tralino lifetime of about T cvt < 2 • As, whereas it does not 
suppress any events with a high lifetime. So, if an event 
contains a neutralino with a long lifetime and which de- 
cays in the detector, the decay photon always has high 
As. 



Here we estimate our sensitivity to neutral, long-lived 
particles which decay to photons in as model-independent 
a manner as possible. We do so by considering a GMSB 
model yJI which we restrict to direct neutralino pair pro- 
duction and decay: pp — > XiXi ~ * 7G7G. We use the 
PYTHIA [3 event generator, with ISAJET [TJ] to generate 
the SUSY masses, and PGS with the parameter file for the 
CDF detector 01 as a simple detector simulation, modified 
for the use of timing information. We accept photons with 
a rapidity \r)\ < 2.1 and a transverse energy Et > 12 GeV 
according to the CDF/EMTiming fiducial region and trig- 
ger 0.01 We first look at the efficiency of the timing system 
with infinite resolution as a function of neutralino mass and 
lifetime for different As restrictions. Then we discuss back- 
ground estimations, take into account a timing resolution of 
1.0 ns and find the sensitivities for the model predictions for 
both the single and diphoton analysis. 

Figure |3 shows the efficiency versus neutralino lifetime for 
a mass of 70 GeV/c 2 [l<| for events with As > ns (pho- 
tons from neutralinos remaining in the detector) and events 
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Table I: The systematic uncertainties, estimated based on 
Refs. for luminosity, acceptance and number of back- 

ground events for use in all analyses in estimating cross section 
limits. 



Factor Syst. 


Uncertainty 


Luminosity 


5% 


Acceptance 


10% 


Number of background events 


30% 



evt 



(ns) 



occurs from a change in the shape of the £z. distribution of 
the neutralinos as shown in Fig.0 For a mass of 80 GeV/c 2 
the maximum moves towards higher ^x. and the distribution 
broadens compared to 40 GeV/c 2 , yielding a greater fraction 
of high-pr neutralinos and hence a loss in efficiency. As the 
mass gets higher the maximum remains the same and the 
distribution narrows, which in turn leads to a gain in effi- 
ciency. Thus, the efficiency is essentially independent of the 
neutralino mass, with slight variations originating from the 
production mechanism, specifically the neutralino momentum 
distribution. 



Figure 4: The efficiency as a function of the event lifetime, 
Tevtj of the neutralino. We distinguish between events in 
which the neutralino remains in the detector, and events with 
photons of medium and large As. The efficiency is 100% 
for prompt decays (a small difference shows up as a binning 
effect) for a photon to be identified, but only a very small 
efficiency for events with low r ov t at large As. At large r ev t 
only few events stay in the detector, however if a neutralino 
is long-lived and stays in the detector, it has large As. We 
note that the true efficiency shape depends on the production 
mechanism i.e. the neutralino pr distribution. 



with a As > 5 ns, separated into single and diphoton events. 
We get essentially the same shapes as in Fig. [1] however the 
overall efficiency is less as we now consider the exponentially 
distributed neutralino lifetime instead of an event lifetime. 
For all four distributions there is an efficiency maximum in 
the lifetime region between 4 and 9 ns. At lower lifetimes 
the probability that the neutralino stays in the detector is 
large enough that the diphoton final state dominates. For 
any As > requirement the efficiency is zero at = ns, 
since all photons have As = 0. At a lifetime of about 3 ns, 
independent of the As cut, single photon events become dom- 
inant. At high lifetimes the efficiency decreases rapidly for 
both analyses as most of the neutralinos leave the detector. 
Hence, in order to have sensitivity in as much lifetime range 
as possible, we consider both 7 + IpT and 77 + IpT analyses. 

In contrast the timing efficiency is essentially constant as a 
function of neutralino mass at a fixed lifetime. Figure[@]shows 
the efficiencies at = 10 ns, where the system has the highest 
efficiency and single photon events dominate. Note that the 
"dip" in the efficiency can be explained by the neutralino pair 
production mechanism: if the decay length is greater than the 
distance to the detector wall, the neutralino will leave. Since 
this is proportional to the ratio of the neutralino's transverse 
momentum to its mass, Sx., (at constant lifetime), the dip 



B. Backgrounds and Sensitivity to Neutralino Pair 
Production 

We now estimate the sensitivity of our system in a quasi- 
model-independent manner using the neutralino pair produc- 
tion introduced in the previous section, but taking into ac- 
count SM backgrounds, more realistic cuts and the timing res- 
olution. We consider separately single photon and diphoton 
events and present our sensitivity as the expected 95% C.L. 
cross section upper limits for either case assuming no signal 
in the data. We also compare our results with the results 
of no EMTiming system available to estimate the effect of 
the EMTiming system over a set of kinematics-only selection 
requirements. Throughout this section we use the relative 
systematic uncertainties for luminosity, acceptance and back- 
ground rates given in Table Q and a As resolution of 1.0 ns. 
The expected cross section limits are calculated following |2fT | 
with the number of events observed "in the data" fluctuat- 
ing around an expected mean background rate according to 
Poisson statistics. The cross section limit is, for a certain 
luminosity, a function of background events and signal accep- 
tance, where both in turn are functions of specified cuts (e.g. 
As and IpT cuts in the 77 + IpT case). By varying the cuts 
we find a signal acceptance and number of background events 
that, after smearing by systematic errors, minimizes the cross 
section limit. 



1. 77 + IpT 

A 77 + IpT analysis is expected to have the best sensitiv- 
ity for low neutralino lifetimes. We follow the analysis in 0] 
(summarized in TableEJ and allow events in which both pho- 
tons have Et > 12 GeV and \rj\ < 1, and study final selection 
requirements on IpT and As. The background for this analy- 
sis consists of QCD events with fake IpT [2l| . We model the 
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Figure 5: The efficiency for events to pass the various As- 
cuts (assuming perfect measurements) as a function of the 
neutralino lifetime at constant mass (m^= 70 GeV/c ), sepa- 
rated into single and diphoton events at As > ns and 5 ns. 
For any As > the efficiency is zero at = ns, since all 
photons would have As = 0. For high r neutralinos have a 
higher probability to leave the detector. For any As one can 
find an efficiency maximum at about 5-10 ns. Single photon 
events are preferred towards higher As requirements and/or 
higher lifetimes, due to increasing probability for a photon to 
leave the detector. Thus, we expect a 77 + |?Tto provide the 
best sensitivity for very low lifetimes, and a 7 + Ipr analysis 
to be best for higher lifetimes. 



Figure 6: The efficiency as a function of the neutralino mass 
at a lifetime = 10 ns for single and diphoton events at As 
> ns and 5 ns (assuming perfect measurements). The ratio 
of single to diphoton events is independent of the neutralino 
mass and is roughly constant as a function of As. One can 
see a soft "dip" in the efficiency curve in a mass range of 
40 GeV/c 2 — 80 GeV/c 2 . This effect is production dependent 
and due to a change in the pr distribution of the neutralinos 
(see Fig. 0. 
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Figure 7: The neutralino ^ distribution for masses 
40 GeV/c 2 , 80 GeV/c 2 and 140 GeV/c 2 for neutralino pair 
production. For a mass of 80 GeV/c 2 the maximum moves 
towards higher ^ and the distribution broadens compared to 
40 GeV/c 2 , yielding a greater fraction of high pr neutralinos 
which either leave the detector or produce low As photons, 
and thus a loss in efficiency. For higher masses the maxi- 
mum remains constant and the distribution narrows so the 
efficiency rises. 



IpT from QCD with a resolution of 10 GeV, i.e. we assume a 
measurement uncertainty of the transverse energy of all par- 
ticles of 10 GeV in each x- and y-direction, as this reproduces 
well the numbers in Q and allows us to extend our search 
region to large values of IpT- Since all photons from QCD are 
promptly produced, we model them with a mean As of ns, 
and a resolution of 1.0 ns. We found that adding the As 
values, Asi2 = Asi + As2, and selecting signal events with 
either large IpT or large Asi2, either of which is not SM-like, 
maximizes the separation of signal and background as shown 
in Fig. [S] The position of the cuts are optimized for each mass 
and lifetime to minimize the 95% C.L. cross section limit, and 
we find that both the Asi2 and IpT cuts are stable at around 
7 ns and 50 GeV for non-zero lifetimes. Without timing 
information we found the optimal IpT cut to also be around 
50 GeV. 



2. 7 + E T + jets 

From efficiency considerations and since the signal does not 
produce jets at the parton level we expect a 7 + IpT + jets 
analysis to yield the best sensitivity for longer neutralino 
lifetimes. We follow the analysis in ^JJ (summarized in 
Table |TTlJ and require the highest Et photon to have 
E T > 25 GeV and \rj\ < 1, a minimum IpT of 25 GeV, and 
no jets or additional photons with Et > 15 GeV, and study 



Table II: The background and baseline selection criteria used 
for the 77 + IpT analysis following Ref. [H.l2ll|. 



Baseline selection requirements: 
E' T X > 12 GeV, E' T 2 > 12 GeV 

< 1, \V J2 \ < 1 
Backgrounds: 

2,577 events / 100 pb' 1 from QCD 

AS12 = As 71 + As 72 , Asi2 = ns with resolution o~As 12 = 
1.41 ns 

IpT- Rayleigh distribution (Square-root of the sum of 2 Gaus- 
sians squared) with a = 10 GeV 

Optimization: 

Accept events where the event has a IpT greater than the 
optimized IpT cut or whose photon has a Asi2 greater than 
the optimized Asi2 cut. 



the final selection requirements on IpT and As. The back- 
ground for this analysis is dominantly QCD, Z7 and cosmic 
ray sources. Since photons from cosmic ray sources hit the 
detector with no correlation between the arrival time and the 
time of collision, we expect them to be randomly distributed 
over time and model this with a flat random distribution in 
As. As in the previous section the As of all other SM sources 
is dominated by the timing resolution of 1.0 ns. The ^rfor the 
backgrounds are modeled according to the shapes in and 
extrapolated to large values using an exponential fit. The ex- 
pected background and signal shapes are shown in Fig.El The 
final cuts on IpT and As are applied to sort out events with a 
large IpT and whose photon has a As within a lower bound and 
an upper bound: Asi < As < As2, to reject photons from 
SM background as well as from cosmic ray sources. We find 
the optimized cuts at around Asi = —2.0 ns, As2 = 2.0 ns 
and IpT = 80 GeV. However, As2 could vary up to 3 ns for 
high lifetimes, IpT up to 120 GeV for higher masses. Without 
timing information we found the optimal IpT cut mostly at 
around 100 GeV. 



3. Results 

Figures IT7T1 and fTTI show the expected 95% C.L. cross sec- 
tion upper limits vs. for = 70 GeV/c 2 and vs. 
for Tx = 20 ns for both analyses for a luminosity of 2 fb _1 . 
Figure IT2*1 shows the ratio of the lower of the two 95% C.L. 
cross section limits with EMTiming system usage and with- 
out in two dimensions. In these plots we see four trends: 1) 
As a function of lifetime the cross section limits rise since the 
probability that the neutralinos decay in the detector goes 
down, 2) at high lifetimes the timing handle is better able to 
separate the signal from the backgrounds and produces better 
limits relative to kinematics alone, 3) as a function of mass 
the cross section decreases as more and more events pass the 
kinematic threshold, and 4) at low masses the timing handle 
is better able to separate the signal from the backgrounds be- 
cause the momentum distribution happens be lower on aver- 
age (see Fig. Eh which shows an example where the kinematics 
are such that there is only small additional acceptance from 
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Figure 8: The distribution of Asi2 and IpT for signal and background in the 77 + IpT analysis. The distributions are (a) from 
direct neutralino pair production, with = 70 GeV/c 2 and = 10 ns, (b) QCD background and (c) direct neutralino pair 
production, with = 130 GeV/c 2 and = 10 ns. The solid and dashed lines show the cuts with and without EMTiming 
system usage respectively that give the smallest 95% C.L. cross section limit. In (c) the mass is so large that there is only small 
additional acceptance from allowing large AS12 events which is why the sensitivity is not improved in this mass region fsee llOfl . 



Table III: The background and baseline selection criteria for 
the 7 + IpT + jets analysis following Ref. [TH |. 



Baseline selection requirements: 
E T > 25 GeV 
W\ < 1.0 
E T > 25 GeV 

No jets or additional photons with E^ 1 > 15 GeV 
Backgrounds: 

12.6 events / 100 pb _i from Z7 ->• uuj, W7, W -> eu, QCD 
and cosmics 

Non-cosmics: As = 0, AS12 = ns with resolution <tas = 
1.0 ns 

Cosmics: 57.2% of total background, flat distribution in As 
IpT distribution according to |l(l |. and extrapolated using an 
exponential 
Optimization: 

Accept events where the event has a IpT greater than the 
optimized cut IpT and whose photon has a As within a range 
of optimized cuts Asi < As < As2. 



allowing large Asi2 events). As expected the 77 + IpT anal- 
ysis yields lower cross sections when the mass or the lifetime 
is low (see Fig. 1121 . The ratio is greatest in this region and 
occurs at a mass around 50 GeV/c 2 and a lifetime of 10-20 ns. 
The 7 + IpT + jets analysis yields lower cross section limits 
for the rest of the considered lifetime range and masses, but 
it is important to note that the course of the separation line 
of the analyses depends on the production momentum distri- 
bution of the neutralino. Unfortunately this analysis cannot 
be applied to searches for long-lived NLSP neutralinos in a 
true GMSB model with the preferred production processes as 
there the neutralinos are produced as part of cascades from 
gaugino pairs which contain jets. Therefore, we do a separate 



analysis for a full GMSB production in the next section. 



IV. SENSITIVITY TO GMSB MODELS 

We next consider the sensitivity to full GMSB production 
where we allow all processes to contribute to the final state 
according to their predicted cross sections. We use the same 
simulation tools as in Section fill Al with the GMSB param- 
eters chosen according to the Snowmass Slope guidelines [l5| 
in the range where the neutralino is the NLSP. Again we con- 
sider a single photon and a diphoton analysis. The 77 + IpT 
analysis methodology is identical to the case of neutralino pair 
production, but the single photon analysis must be modified 
to allow jets as here the neutralinos are part of cascades from 
gauginos which produce additional particles which, in general, 
could be identified jets. We thus use a 7 + IpT + jets analysis. 



A. 7 + Et + jets 

A 7 + Ipr + jets analysis should be most sensitive to neu- 
tralinos with long lifetime which are produced in association 
with other particles in the final state such as from gaugino 
pair pro duct ion. We follow the analysis in 0] (summarized 
111 Table ITv) and require events with the primary (highest 
E T ) photon to have E T > 25 GeV and = [0, 1.1] or [1.5, 
2.0], E T > 25 GeV, > 2 jets of E T > 20 GeV and < 2.0, 
and study the final selection requirements on IpT and As. 
The backgrounds are dominated by QCD and W+jets |22T |. 
The expected IpT of the background is modeled according to 
Ref. and since the backgrounds are from SM we take a 
mean As = ns with a resolution of 1.0 ns. The signal and 
background shapes are shown in Fig. ED We find that the op- 
timal final selection requirements accept events in which the 
event has a large IpT or a large As. Again, both cuts are op- 
timized to minimize the 95% C.L. cross section limit for each 
mass and lifetime case. For without-timing usage we find the 
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Figure 9: The distribution of As and ~$t for signal and 
background in the 7 + JpT + jets analysis. The distri- 
butions are (a) from direct neutralino pair production, with 
= 70 GeV/c 2 and = 10 ns, and (b) from SM back- 
ground. The solid and dashed lines show the cuts with and 
without EMTiming system usage respectively that give the 
smallest 95% C.L. cross section limit. 
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Figure 10: The expected 95% C.L. cross sections lim- 
its on neutralino pair production as a function of for 
nix = 70 GeV/c 2 (top) and as a function of for = 20 ns 
(bottom) in the 77 + Ifir analysis for 2 fb -1 luminosity both 
with and without a timing system for comparison. As ex- 
pected at = ns the cross sections merge as the timing 
system has no effect; for higher the sensitivity goes down 
as more photons leave the detector, but the difference of the 
limits increases as As gets larger for the signal and timing 
becomes more helpful. The limits get better as the mass goes 
up since more of the events pass the kinematic requirements, 
however the timing system only provides real additional sen- 
sitivity at the lowest masses where the neutralino momentum 
distribution is softer. 



IpT cut to be around 50 GeV for masses around 70 GeV/c 2 , 
varying up to 110 GeV for masses around 150 GeV/c 2 . For 
with-timing usage we find only a As cut around 4 ns which is 
stable for all masses and lifetimes, and no I$t cut other than 
the baseline IpT > 25 GeV (except for = ns where the 
diphoton case has the best sensitivity). While it is outside 
of our ability to predict, one might find further optimization 
is possible by further lowering the baseline selection require- 
ments. 



B. Results 

Figures ITU and ITol show the expected 95% C.L. cross sec- 
tion upper limits vs. for = 70 GeV/c 2 and vs. 
for = 20 ns for both analyses for a luminosity of 2 fb _1 . 
Figure EJ shows the ratio of the lowest 95% C.L. cross sec- 
tion limits with EMTiming system usage and without in two 
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Figure 11: The expected 95% C.L. cross section limits on neu- 
tralino pair production as a function of for = 70 GeV/c 2 
(top) and as a function of for = 20 ns (bottom) in the 
7 + IpT + jets analysis for 2 fb _1 both with and without a 
timing system. As in Fig. 1101 for higher the sensitivity goes 
down as more photons leave the detector, but the difference 
of the limits increases as As gets larger for the signal and 
timing becomes more helpful. The curves do not merge at 
ns lifetime since cosmic ray backgrounds always contribute 
at high As and the timing system always has some effect on 
the cross section limit. The rise from 10 ns to ns originates 
in an increasing probability towards zero lifetime for two pho- 
tons to remain in the detector, yielding lower efficiency for a 
single photon analysis. The limits get better as the mass goes 
up since more of the events pass the kinematic requirements. 



dimensions. We see the same general trends as in neutralino 
pair production as the signal shapes are similar in both anal- 
yses. 

A comparison of the cross section limits with the produc- 
tion cross sections in the GMSB model at a luminosity of 
2 fb _1 gives the mass vs. lifetime exclusion regions shown 
in Fig. 1171 As expected, timing has the biggest effect at low 
masses and high lifetimes. We have also indicated the ex- 
clusion regions from LEP II from both direct and indirect 
searches [13| . and the line = 1 keV/c 2 , below which 
is the theoretically preferred region from cosmological con- 
straints [T3 |. LEP effectively excludes all neutralino masses 
under 80 GeV/c 2 up to high lifetimes, with a small exten- 
sion to 100 GeV/c 2 for lifetimes below 20 ns. For 2 fb -1 , in 
run II, the Tevatron should significantly extend the sensitiv- 
ity at large mass and lifetimes, covering most of the lifetimes 
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Figure 12: This plot combines the 77 + Ipr and 
7 + JpT + jets analysis results for neutralino pair produc- 
tion for 2 fb -1 of data and is a 2-dimensional visualization of 
Figs. llO| and llll The contours of constant cross section limit 
are shown as the solid lines, and the separation line between 
the regions where the two different analyses provide the best 
sensitivity is given by the dotted line; the 7 + $rf jets 
analysis shows better cross section limits than a 77 + Jfir 
analysis in the mass and lifetime range above the dashed line. 
The shaded regions delineate the contours of the ratio of the 
95% C.L. cross section limits between with and without EM- 
Timing information. The ratio is greatest for a low neutralino 
mass and a lifetime of 10-20 ns, and lowest for a high mass 
and low lifetime. 



Table IV: The background and baseline selection criteria used 
for the 7 + IpT + jets analysis following Refs. [Oi US- 



Baseline selection requirements: 
E% > 20 GeV 

|?7 7 | < 1.1 and 1.5 < |?7 7 | < 2.0 
IJjt > 25 GeV 

>2 jets with E J T et > 20 GeV and \rf et \ < 2.0 

Backgrounds: 

320 e vents / 100 pb" 1 from QCD and W+jets 

AS12 = ns with resolution <ta s = 1.0 ns 

Ipr distribution from (Tll |. extrapolated to large IpT 

Optimization: 

Accept events where the event has a Ipr greater than the 
optimized Ipr cut or whose photon has a As greater than the 
optimized As cut. 



for uiq < 1 keV/c up to a mass of around 150 GeV/c . The 
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Figure 13: The distribution of As and Ipr for signal and back- 
ground in the 7 + IpT + jets analysis. The distributions are 
(a) from full GMSB production, with = 70 GeV/c 2 and 
Tx = 10 ns, and (b) from SM background. The solid and 
dashed lines show the cuts with and without EMTiming sys- 
tem usage respectively that give the smallest 95% C.L. cross 
section limit. 



mass exclusion limit at 168 GeV for = ns is comparable 
to the limit presented in the D0 study of displaced photons 
in Ref. [zsTj . but for large lifetimes this result significantly ex- 
tends the reach. 
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Figure 14: The expected 95% C.L. cross section limits as a 
function of for = 70 GeV/c 2 (top) and as a function of 
rrix for = 20 ns (bottom) at 2 fb" 1 for with and without 
EMTiming for full GMSB production in a 77 + IpT analysis. 
The results are similar to those in Fig. 1101 



Factors that might change the cross section 
limit 



While we have taken the best available nominal values from 
the references for both the contamination of cosmic ray back- 
ground events and the timing resolution in the 77 + IpT and 
a 7 + IpT + jets analyses, the limits are sensitive to mis- 
estimations of these values. For simplicity, rather than in- 
clude them as a systematic error we estimate the variation 
of our results for these effects on our cross section limits for 
a neutralino mass of 110 GeV/c 2 and a lifetime of 40 ns be- 
yond the LEP II exclusion region. Figure Hsl shows the cross 
section limit as a function of the fraction of events which are 
from cosmics in the background sample. Using the same anal- 
ysis style as in the 7 + IpT + jets case we find cuts around 
Asi = 3.0 ns and IpT = 55 GeV, with As2 varying from in- 
finity down to 7 ns with a rising fraction of cosmics. The 
limits rise approximately linearly as a function of the fraction 
of events which are from cosmics. An upper bound on the 
fraction of cosmics of 10% would reduce the limits by about 
a factor of four; a more reasonable estimate is probably 1-5% 
which would raise the limts by a factor of 2-3. The limits are 
potentially more sensitive to the resolution. Figure ED shows 
how the limits change as a function of the system resolution 
for the same mass and lifetime, for the original 7 + IpT + jets 
analysis. While there is a mass/lifetime dependent resolution 
threshold from where the limit can change drastically and 
approach the cross section without EMTiming, the limits are 
fairly stable (with the same factor of 2) for resolutions within 
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Figure 15: The expected 95% C.L. cross section limits as 
a function of for = 70 GeV/c 2 (top) and as a func- 
tion of rnj for = 20 ns (bottom) at 2 fb~ for with and 
without EMTiming System for full GMSB production in a 
7 + JfiT + jets analysis. For all but the lowest lifetimes the tim- 
ing information significantly improves the cross section limits. 
Note that here the cross sections merge at zero lifetime since 
we have neglected cosmics in this analysis following [Tll.l22T|. 



20% of the nominal 1 ns resolution, and there is good reason to 
believe that the resolution will be better than advertised Q. 



Figure 16: This plot combines the 77 + ]pr and 7 + $t + jets 
analysis results for a full GMSB model simulation for 2 fb _1 
of data and is a 2-dimensional visualization of Figs. 1141 and 
fTTH The contours of constant cross section limit are shown 
as the solid lines, and the separation line between the regions 
where the two different analyses provide the best sensitivity 
is given by the dotted line; the 7 + I£t + jets analysis shows 
better cross section limits than a 77 + I£t analysis in the 
mass and lifetime range above the dashed line. The shaded 
regions delineate the contours of constant ratio of the 95% 
C.L.- cross-section limits between with and without EMTim- 
ing information. The EMTiming system has its most effective 
region at high lifetime and is least effective for high masses, 
where the kinematics give the best separation. 



V. CONCLUSION 

We have studied the prospects of using the new EMTiming 
system at CDF in the search for neutral, long-lived particles 
which decay to photons, as one can find in supersymmetric 
models. We find that the kinematic requirements and the EM- 
Timing system provide excellent rejection against SM back- 
grounds in complementary fashion. As the mass increases the 
kinematics are more important and the sensitivity gets bet- 
ter. For a given mass, as the lifetime increases more and more 
of the neutralinos leave the detector and the overall sensitiv- 
ity goes down, but the EMTiming system provides additional 
rejection power and allows for significant exclusions even at 
large lifetimes. While the region where EMTiming produces 
the most additional rejection is already excluded by LEP II, 
the additional handle it provides should allow the Tevatron in 
run II to produce the world's most stringent limit at masses 
above 80 GeV/c 2 at high lifetimes and has the potential to 
cover the entire region for mg < 1 keV/c 2 up to a neutralino 
mass of around 150 GeV/c 2 . 
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Appendix A: PHOTON POINTING 

As shown in Fig. El LEP II has already excluded the low 
neutralino mass region using a photon pointing method [lI |. 
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Figure 17: The expected 95% C.L. exclusion regions as a func- 
tion of neutralino lifetime and mass for full GMSB production 
at 2 fb -1 luminosity for the the 77 + Ipr and a 7 + $t + jets 
analysis separately. In plot (a) the region below the dashed 
line is the exclusion region from kinematics alone, i.e., where 
no timing information is used. Plot (b) shows the full ex- 
clusion region from the overlap of the two analyses and com- 
pares the results to the direct and indirect search limits from 
LEP II [lI | and the theoretically favored region from cosmo- 
logical constraints (mg < 1 keV/c 2 ) 01- The Tevatron in 
run II should be able to significantly extend the LEP II limits 
and provide sensitivity in the favored region for all masses 
below about 150 GeV/c 2 . 



Figure 18: The 95% C.L. cross section limit on full GMSB 
production at = 110 GeV/c 2 and = 40 ns as a function 
of the fraction of the background from cosmic ray sources 
for a 7 + TpT + jets analysis. The cross section limits rise 
approximately linearly as a function of the fraction and 10% 
provides an outer bound on this fraction. A more reasonable 
fraction is probably 1-5% which roughly doubles the limit. 



In this section we compare the EMTiming system to a poten- 
tial photon pointing ability at CDF. A non-zero lifetime can 
result in a macroscopic decay length and impact parameter, 
where the impact parameter of the photon is basically the 
closest distance of the trajectory to the collision point. While 
CDF has never used its calorimeter for a pointing measure- 
ment it is possible to use the central EM strip/wire gas cham- 
ber (CES) and the central pre-radiator gas chamber (CPR) 
at CDF to measure two points along the photon trajectories 
that determine the direction of the photon, and trace it back 
to yield the impact parameter j^J- Since the CPR has no 
^-measurement ability this allows only a measurement of the 
radial component of the impact parameter with an estimated 
resolution of 10 cm (see Table lvjl . One of the primary reasons 
this has not been used is the conversion, i.e. measurement, 
probability, of ~65%, with an angular dependence obtained 
with: 

__md 

Pc = 1 - e , 

where 7V ra d = 1.072 is the number of radiation lengths before 
the CPR and 6 is the angle with respect to the beamline. 
To estimate the sensitivity with a pointing method we con- 
sider a 7 + IpT + jets analysis. Figure [20l shows the distribu- 
tion of the signal events vs. impact parameter and As in a 
7 + IpT + jets analysis taking into account the measurement 
probability; there are roughly as many events in the region of 
low impact parameter and high As as there are at high im- 
pact parameter and low As. Hence either method should have 
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Figure 19: The 95% CL. cross section limits vs. EMTim- 
ing system resolution for = 110 GeV/c 2 and r-j = 40 ns 
at 2 fb -1 luminosity in the 7 + Ifir + jets analysis. As ex- 
pected, for large resolution the cross section with EMTiming 
approaches the cross section without EMTiming. A system 
resolution of 1.0 ns improves the cross section limit by a factor 
of about 20, but this varies as a function of mass, lifetime and 
depends on the analysis. It is reasonable to assume that the 
resolution will be within 20% of the nominal 1.0 ns presented 
here. 



roughly the same effect on the exclusion region, as confirmed 
by Fig. |^ which shows the expected exclusion region in the 
mass-lifetime plane. While timing is better than pointing by 
itself, if pointing turns out to be feasible, a combination of 
the two would improve the sensitivity. 



Considered separately, a second advantage of timing is that 
it "filters" manifestly long lifetime events, whereas the impact 
parameter allows also short lifetime-high momentum events, 
which might come from SM. Another possible advantage of 
the combination is that in the event of an excess, we could 
draw more information about the individual events, for in- 
stance determine the direction of the photon. With the x- 
y-direction being fixed by the CPR/CES measurement, we 
can use the timing system to measure the ^-direction. Or if 
the pointing provided the photon direction in z and x-y, one 
could possibly determine the position of the vertex and thus 
the decay time. However, with the current 1.0 ns resolution 
the photon vertex position resolution would be roughly 50 cm, 
if we assume the neutralino boost to be ~ 1.0. 



Table V: Photon pointing parameters for the CDF detector 
[24|]. With this combination we estimate an impact parameter 
resolution of 10 cm in the radial direction. 
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Figure 20: The relationship between As and the impact pa- 
rameter 6 of a photon from x? — > 7G decays in a GMSB model 
with = 70 GeV/c 2 and = 10 ns. We show the selection 
requirements that give us the smallest 95% C.L. cross section 
limit in a 7 + IpT + jets analysis. The photons without impact 
parameter measurement are assigned a b = —0.1 m. Due to 
the low cut on the impact parameter there are about as many 
events in the low-As high-6 as in the high-As low-b quarter, 
leading to a similar efficiency for a pure b-cut compared to 
a pure As cut. The combined restriction leads to improved 
signal sensitivity. 



14 



-£300 

CD 
E 
0) 

= 250 

o ,_ 



200 



150 



100 



50 



' I 1 1 1 I 1 1 1 I 1 

Predicted exclusion region with £ 
y+^ T +jets analysis with timing 

Predicted exclusion region 
with pointing 

Predicted exclusion region 
with pointing and timing 



2 fb luminosity 



± 



60 80 100 120 _ Q 140 160 

X 1 mass (GeV/c ) 



Figure 21: A comparison of the expected exclusion regions 
as a function of neutralino mass and lifetime for the GMSB 
model at 2 fb -1 luminosity for a 7 + $t + jets analysis with 
photon pointing and timing. While timing generally yields 
a higher sensitivity than pointing, both methods would, if 
we combined them, extend the exclusion region further than 
either of them alone. 
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